Another Chirality Quantum Phase Transition in the Atom-Cavity Interaction 
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This paper is devoted to a study of a kind of chirality appearing in a model in cavity or circuit 
quantum electrodynamics. The model especially describes an artificial atom fully coupled to a 1- 
mode photon in a cavity when the interaction is in ultra-strong coupling regime. This chirality is 
played in the fully coupled (FC) Hamiltonian by two modified Jaynes-Cummings (JC) Hamiltonians. 
The modified JC Hamiltonian as well as the standard one has the superradiant ground state energy 
for sufficiently large coupling strength. The superradiance is caused by a phase transition found 
by Preparata, called the Dicke-Preparata (DP) superradiance. Using the one common Hamiltonian 
formalism with the chirality, we theoretically prove the followings: (i) when the coupling strength 
makes the weak or standard strong coupling regime, the ground state energy of the FC Hamiltonian 
can be well approximated with that of the standard and modified JC Hamiltonians since the chirality 
does not work in the regime; (ii) once the coupling strength plunges into the ultra-strong coupling 
regime, the effect of the chirality is turned on and works in the ground state energy of the FC 
Hamiltonian. That is, there is a kind of a phase transition concerning the chirality associated with 
the DP superradiance between the strong coupling regime and the ultra-strong coupling one. 

PACS numbers: 42.50.-p,03.65.Ca,03.65.Db,85.25.-j 



I. INTRODUCTION 

The interaction between an atom and the light in na- 
ture follows the quantum electrodynamics (QED), and 
its coupling strength is governed by the fine structure 
constant. Thus, the perturbation theory is valid over the 
standard QED world. On the other hand, cavity QED 
supplies stronger interaction than the standard QED 
does 

QH3. Namely, it can realize the strong coupling 
regime for the atom and the light in a cavity. Such a 
strong interaction is usually prepared with a two-level 
atom (or qubit) in a mirror cavity (i.e., a mirror res- 
onator) . Several solid-state analogues of strong coupling 
regime of cavity QED have been foreseen in supercon- 
ducting systems (3-fl2|. To sum up, we respectively re- 
place the two-level atom, the laser, and the mirror res- 
onator by an artificial atom, a microwave, and a mi- 
crowave resonator on a superconducting circuit. This 
replaced cavity QED is called circuit QED. Actually, it 
has been demonstrated experimentally [l3l - l20| . It is re- 
markable that circuit QED is capable of realizing the 
so-called ultra-strong coupling regime for the interaction 
between the atom and the light [2ll - |24l ]. Namely, cir- 
cuit QED has given us the technology which implements 
much stronger atom-light interaction than the standard 
QED does. Although the Jaynes-Cummings (JC) model, 
namely the rotating wave approximation (RWA), is use- 
ful to analyze several physical properties in the weak 
and the standard strong coupling regimes, it has been 
reported that JC model is no longer valid over the ultra- 
strong coupling regime and that the effect coming from 
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the counter-rotating terms is important [2J-|32[. There- 
fore, there is an urgent need to make an approximate 
Hamiltonian formalism [33l [34[ which works as a non- 
perturbative theory for the ultra-strong coupling regime 
in circuit QED. Then, the points are that the theory 
should explain the effect of the counter-rotating terms 
and that the approximate Hamiltonians should be solv- 
able. 

This paper deals with a model in the theory of cir- 
cuit QED, which describes the energy of a 2-level arti- 
ficial atom fully coupled to a 1-mode photon in a cav- 
ity. We call its total Hamiltonian the fully coupled (FC) 
Hamiltonian in this paper, though it is also called the 
Rabi Hamiltonian. In Ref.[35| Irish has proposed the 
generalized rotating wave approximation (GRWA). The 
GRWA comes up with the same expression of eigenstate 
energies as given by Feranchuk et al. [3fjj . In numeri- 
cal analysis, the GRWA is valid over even ultra-strong 
coupling regime. This expression suggests that a modi- 
fied JC model should be useful for the eigenstate-encrgy 
problem for the FC model. Meanwhile, the modified JC 
Hamiltonian as well as the standard one has the superra- 
diant ground state caused by the Dicke-type energy level 
crossings (37rl42l |. These crossings take place in the case 
where the coupling strength is sufficiently large. This 
type of superradiance is the Dicke-Preparata's [HI, |H| . 
That is, it is a kind of phase transition so that the 2-level 
system coupled to a photon-field switches from the ini- 
tial ground state to a new non-perturbative ground state 
[HI, [3^ . We decompose the FC Hamiltonian into a mod- 
ified JC Hamiltonian and its chiral-counter Hamiltonian 
in the same way as in Ref.[38|. The former acts in the 
standard state space and the latter is among the modified 
JC Hamiltonians in the chiral state space. But, since the 
chiral-counter Hamiltonian controls the counter-rotating 
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terms in the standard state space, it breaks the proper- 
ties special to the modified JC Hamiltonian. We call the 
pair of the former and latter the chiral pair. Using the 
mathematical structure of the chiral pair, in the same 
way that Bermudez et al. have studied a quantum phase 
transition in the Dirac oscillator from the point of the 
view of the chirality [ill H(| , we will show that the FC 
Hamiltonian has the phase transition concerning the chi- 
rality associated with the Dicke-Preparata (DP) super- 
radiance when the atom-cavity system plunges into the 
ultra-strong coupling regime from the standard strong 
coupling one. We also show how the effect of the chiral- 
ity increases as the coupling strength grows large. In fact, 
the phase transition associated with the DP superradi- 
ance appears also in the standard state space when the 
coupling strength approaches the deep-strong coupling 
regime, but we realize that it takes place faster in the chi- 
ral state space than it does in the standard state space. In 
addition, by using our decomposition and estimating the 
non-commutativity between the two component Hamil- 
tonians of the chiral pair, we will prove that the chirality 
associated with the DP superradiance does not work in 
the weak and standard strong coupling regimes, and thus, 
the ground state energy of the FC Hamiltonian is well 
approximated by that of the (modified) JC Hamiltonian. 
Some effective Hamiltonians are individually studied in 
Refs.|33l.l34j to handle the weak, strong, and ultra-strong 
coupling regimes. We add another observation to them 
through our decomposition to handle all regimes in the 
one common Hamiltonian formalism. 

Our paper is constructed as follows. In Sec|TT]we re- 
call some facts on the ground state energy of the FC 
Hamiltonian, and formulate our problems. In Sec lIIH 
parameterizing the standard JC Hamiltonian, we define 
a modified JC Hamiltonian and its chiral-countcr Hamil- 
tonian to handle the counter-rotating terms. In Sec lIVl 
we recall the Dickc-typc energy level crossing and the 
superradiant ground state energy for the modified JC 
Hamiltonian. In Secs|V] and I VII we give answers to our 
problems formulated in Seclrfl 



II. HAMILTONIAN OF FULLY COUPLED 
MODEL 

In this section, we recall a mathematically exact ex- 
pression of the ground state energy of the fully coupled 
(FC) Hamiltonian, and formulate our problems. Hamil- 
tonians appearing in cavity QED and circuit QED fas- 
cinate some mathematicians. They belong to a class of 
matrix- valued Schrbdinger operators |47j . 



A. Fully Coupled Hamiltonian 

We consider the Hamiltonian describing the energy of 
the 2-level atom fully coupled to a 1-mode photon. We 
denote the creation (resp. annihilation) operator for the 



1-mode photon by cl< (resp. a). We use the standard 
notation for the the Pauli matrices as a x = ( \ \ ) , <J y = 

In addition, we denote the 
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2x2 identity matrix by ctq, i.e., ctq = ( J \ 



Our FC Hamiltonian Hpn reads 



tlFC '■= ~^~ a z 



a>a+- 



hga x (c 



This is the 1-mode photon version of the so-called spin- 
boson model [38l440l l48l - [5l| . Here the constants oj a and 
lo c are respectively the atom transition frequency and 
the cavity resonance frequency, and the parameter g > 
stands for the atom-photon coupling constant. Follow- 
ing Ref . j3l| . in this paper, we temporarily define the 
ultra-strong coupling regime by the region where the di- 
mensionless coupling constant g/w c satisfying the condi- 
tion 0.1 ^ g/u> c . Thus, we regard the region with the 
condition g/w c ^ 0.1 as the weak and strong coupling 
regimes. In addition, we call the region with the condi- 
tion 1 ^ g/w c the deep-strong coupling regime. 

We denote by E-pc the ground state energy of the 
FC Hamiltonian H-pc- Since the FC model is the 1- 
mode photon version of the spin-boson model, we can 
apply several results on the spin-boson Hamiltonian to 
our FC Hamiltonian. In Ref.[5(| we gave a strict ex- 
pression of the ground state energy of the spin-boson 
model using the Bloch formula, the Duhamel formula, 
and the parity conservation. We apply the expression to 
the ground state energy Epc now. To see the expres- 
sion, we introduce some notations. First, we note that 
the FC Hamiltonian Hpc has the following parity con- 
servation: [Hpc, <?z{— 1)^1 = 0. Using this parity con- 
servation, there are many attempts to develop the accu- 
racy of the numerical anal ysis for the eigenvalue problem 
for the FC Hamiltonian |35l . l36l . l52l [53j . The expres- 
sion of the ground state energy Epc has the two parts 
with even photon-numbers and the odd photon-numbers 
respectively. This division comes from the parity con- 
servation. The method for the expressions in Ref. [50| 
reminds us of the computations of the transition ampli- 
tude of the so-called instanton gas by the Euclidean path 
integral fH3 |. 

We define functions I even (f3) and IoddW) of a variable 
/3>0by 



/even(/3) - = ^ + J2 



Jo 



n ■ ■ ■ 



df3 2t 



x e -2(g 2 /Wc)( 2 G/3i,- -,/3 2 «+2^) 



W/3) :=/3^ e - 2 fe 2 /^) 



j2L /lj -.21+1 rP rPi rPn 

+E(^) M dfc--- 4W 

g=1 v 1 ' Jo Jo Jo 

x e -2(g 2 /" c 2 )(2G 01 ,.. ^ 2l +2F^...^ 2t+1 +{U+l)) 



for the sequences {Gp lr .. y /3 2e }^ =l and {Fp lt 



-> oo 
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given by 
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Pi,-, 02 



E e ~ 
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{P2p-1—P2p)w c 
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P,g=l;p<« 



and 



E 

p=l 



-Pzp-lUo 



0<o. 



Then, Theorem 1.3(i) of Ref.[50| says that the ground 
state energy Epc is expressed as 



Efc = - — lim % In {I e 

2 0J C P^roo jj 



1 

(2.1) 

for arbitrary coupling constant g provided that 1/2 < 
Lo c /u> a . We note that this expression was primarily proved 
without the assumption 1/2 < w c /u; a , but under the 
assumption that the vector ip+ defined below overlaps 
with the ground state t/'-Efc °f the FC Hamiltonian i.e., 
(V , BfcIV'+) 7^ 0. Here the vector ip + is given by ip+ '■- 



cos {gP} 1 1, 0) + i sin {gP} | |, 0) and P = 



7.t _ 



The vectors 1 |, 0) and | f, 0) are respectively defined by 
| |,0):=| |) <g> |0) and | f,0):=| t) ® |0) for spin states, 
ji) = (J) and |t) = (J), and the Fock vacuum |0) of 
the photon-field. Theorem 1.3 (ii) of Ref.[5(| says that 
this condition holds if 1 — e~ 2g "/" c < 2uj c /uj a . Then, it is 
always satisfied for all coupling constant g if 1/2 < oj c /ui a . 
Thus, wc assume the condition 



1 U! r 

- < — 

2 oj h 



< 1 



(2.2) 



throughout this paper, where the latter inequality makes 
our arguments after Sec IIIII simple. 

Eq. (|2.1[) is strict, but the expression is very com- 
plicated because those of the functions I cvcn (f3) and 
lodd(P) are so. Thus, we can make it simpler 
with a constant. We modify the functions /even(/9) 
and iodd(/3) by replacing the constants Gp u ...^ 2t and 
Fp lt ... t p 2e+1 in them with IG and G/2, respectively, for 
an arbitrary constant G in the closed interval [—1,0]: 



4ven(/3) := COsh 



(/W2)e 



(/^ a /2)e- 2 s 2 ( G+1 )/-c 



-2g 2 (G+l)/w 



and 1^(3) := 
Then, Theorem 1.5 of 



sinh 

Ref.[]50j says that for every coupling constant g we can 
uniquely determine a constant G(g) in the closed interval 
[—1,0] so that the ground state energy Epc turns out to 



be a simple expression: 

hue hg 2 



Z D c p— >oo p L J 



hide 



-exp 



-2^ (G(g) + 1) 



(2.3) 



The constant G(g) is determined as a solution of the 
equation: 



lim 



/even(/3) - W/3)) 



VP 



r G(g) 



(/3)-4 G c ( d S; (/3) 



1. 



In fact, we have the limit G(g) — > as g — > oo (See 
Sec for its proof). We note that the constant hG{g) in 
Eq. (|2.3p plays a role similar to the classical action asso- 
ciated with a single- instanton solution (see Eq.(3.36) of 
[53 |) in the expression of the transition amplitude. Thus, 

/e™n(/3iand I^fiP) correspond to individual Eq.(3.41) 
of Ref.[H3. 

By taking —1 and as the constant G(g) in Eq. 
we can derive the roughest estimates as in Figflja) 



eiow(g) 



~2~ ~ ~ 

< Epc 

ttiVc _ Hg^ _ _ 2g 2 /u> 2 

- 2 ' uj c 2 6 



e upp (g). (2.4) 



We note that -Efc < huj c /2 — hg 2 /u} c +huj a /(uj a/ /2~g 2 /L0 c ) 
provided SuJaUr. < 4g 2 following the inequalities on p. 161 
of Ref . [55l. IBa] . but we use the upper bound e upp (g) in this 
paper. We have the inequalities < e upp (g) — ei ow (g) < 
hw & /2 and the limit as in FigQJb): lim g _ i . 00 (e upp (g) — 
eiow(g))/fi^a = 1/2. Namely, the difference between the 
two bounds, e upp (g) and ei ow (g), is the zero-point energy 
(i.e., the vacuum fluctuation) at most. 

Let us point out two remarks here. The first one is 
that we have the asymptotic behaviors [56j | : 



Eyc ~ -h(u a - lu c )/2 

if the coupling constant g is so small, 
hm c _ Hg 2 

2 U! c 

if the coupling constant g is so large, 



Efc 



(2.5) 



E 



FC 



'upp 



(g) 



if the ratio 



is so large 



(For the detailed proofs of them, see Sec[5}. 

The other remark is concerned with the second asymp- 
totic behavior of Eqs. (|2.5| ). This limit is due to the 
fact that each energy of the Hamiltonian i/pc is dom- 
inated by that of the Hamiltonian H AA when the ratio 
g/uj a is so large. Here the Hamiltonian H AA is given 
by H AA := fruj c (a^a + |) + hga x {a i + a). This Hamil- 
tonian H AA is unitarily equivalent to the Hamiltonian 
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FIG. 1. Set u> — u) a — u c . (a) e upp (g)/7iu;(red solid line), 
ei ow (g)/?iw(grcen solid lino), numerically-calculated ground state en- 
ergy Efc (blue dotted line); (b) (e upp (g) — ei ow (g)) / hw (black solid line) . 

t\juj c {a) a+^) + Hg<j z (a 1i + a), and thus, unitarily equivalent 

to the Hamiltonian cro {fio;c(a^a+ \) — (^g 2 / w c)}- The 
existence of the identity matrix <tq causes each asymp- 
totically degenerate state for infinitely large coupling 
strength (i.e., g/w a — > oo) as in FigJ5] which has the nu- 
merically computed energy levels of the FC Hamiltonian. 
In particular, Fig(2] says that the first excited state en- 
ergy asymptotically approaches to the ground state one. 
But, we can realize that there is no strict coincidence be- 
tween the two energies because the ground state of the 
FC model is unique for every coupling constant g, apply- 
ing Theorem 3 of Rcf. on the spin-boson model to our 
FC Hamiltonian. We can show the asymptotical degen- 
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strength g/co of interaction 



FIG. 2. Energy Levels of H^c f° r — w a — w c . Each color indicates 
the nth level of the energy for n — 0, 1, 2, • • • from the bottom, where 
the 0th level energy means the ground state energy. We will make a 
remark on the accuracy of the numerical ana lysis for energy levels of 
the FC Hamiltonian at the tail end of Scc lVIl 

eracy of the ground state. That is, for the first excited 



state energy E^ of the FC Hamiltonian, we have the 
asymptotic behavior: 

Sec - E FC = hg x o(^/g 2 ) as c^/g -> (2.6) 

(See SecfA]for its proof). The degenerate ground state 
in circuit QED has been studied for a model of a chain 
of N superconducting Josephson atoms in Ref . (57| . 



B. Formulation of Our Problem 

As far as the coupling constant g makes the weak or 
the standard strong coupling regime of circuit QED, the 
Hamiltonian H F q is well approximated by that of the JC 
model: iff c := ^a z +hu c (a^a + \)+hg (oV_ + a+a), 
where ct_ := (a x — ia y )/2 = (?q) ano - a + :== ( a x + 
i(j y )/2 = ( o o ) ■ This model is the 1-mode photon version 
of the Wigner-Weisskopf (WW) model called in mathe- 
matics and mathematical physics. That is, the JC (resp. 
WW) Hamiltonian is obtained by applying the RWA 
to the FC (resp. spin-boson) one. The WW Hamil- 
tonian has a symmetry for a slightly modified photon- 
number operat or ( see, for instance, Eq.(6.2) of Ref.[49| 
or Eq.(2.21) of |3g])- Mathematically, this symmetry re- 
duces the eigenvalue problem for the JC Hamiltonian 
to the problem of solving a quadratic equation in each 
cigenspace of the modified photon number operator, so 
that the JC Hamiltonian becomes solvable. But the FC 
model docs not have the symmetry. Thus, the eigenvalue 
problem for the FC model has not been solved exactly 
yet, or rather, even whether it is a solvable model has 
not been settled, though Reik and Doucha conjectured 
[EH that the solution of the eigenvalue problem would be 
exactly represented in Bargmann's Hilbcrt space (59j . 

We denote the ground state energy of the JC Hamilto- 
nian by Ej C . As shown in Fig|3la) (where we set w a and 
lo c as ui a = lj c = uj), the ground state energy Ej C can 
supply a good approximation for the upper bound e upp (g) 
until about g/uj w 0.06. Here, denote the sum of counter- 
rotating terms by W(g) := fefaVj. + cr_a). Then, em- 
ploying the idea used in Refs.p2l. |60| . we consider the 
Hamiltonian Hpcil) with a dimensionless parameter 7 
now: H FC (j) ■= #f c + lW(g) for < 7 < 1. We note 
that i?pc(0) = Hf c and H FC {1) = H FC . When the 
parameter 7 is so small that the counter-rotating terms 
7W(g) is perturbative for the JC Hamiltonian Hf c , the 
Hamiltonian Hfc (7) nas the same energy-spectral prop- 
erties as the JC Hamiltonian does [4l|, |42|. We note 
that the introduction of such a parameter was done in 
Ref. [4t| . too. On the other hand, when the parameter 
7 is not so small, we cannot ignore the effect caused by 
the counter-rotating terms jW(g) as well as the impor- 
tance of the effect is claimed in Refe.0 H, For 
instance, let us take 1 as 7, then we have our desired 
Hamiltonian: ffpc(l) = ^fc- As in Figs|31 the approx- 
imation of the ground state energy E F q by the ground 
state energy £j C does not function well any longer after 
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a coupling regime. Therefore, we need another device to 
handle such a coupling regime. Meanwhile, FigJ3{b) says 
that the JC model has the phase transition at g = ui. 
That is, the superradiant ground state energy appears 
as follows: Ej C = (a constant) if < g/oj < 1; 
£j C = hcu — hg (a first-degree polynomial function of g) 
if 1 < g/ w < 1 + v2- We note here that the ground state 
is dressed with no photon in the case < g/w < 1, but it 
is dressed with a photon in the case l<g/w<l + \/2- 
We will precisely recall in Scc lIVI how the superradiant 
ground state energy appears. 




0.02 0.04 0.06 0.08 0.1 0.5 1 1.5 



FIG. 3. Sot uj — cj a — uj c . e upp (g)/?iu;(red solid line), 
e] ow (g)/?iu;(grccn solid line), and E^ G /hco (dark-blue dotted line). 

Introducing another dimensionless parameter e with 
< e < 1, we parameterize the two frequencies u; a and 
lo c , and then, denote the parameterized frequencies by 
w a (e) and w c (e). Replacing the original frequencies in 
the JC Hamiltonian iJjc with the parameterized ones, 
we obtain a family of the modified JC Hamiltonians 

+ hg (a T cj_ + (T+a) . (2.7) 

In Sec lIIIi we will show that there is a parameteriza- 
tion so that the Hamiltonian H-pc is divided into the two 
parts: 

H FC = Hf c (e) + sa x H^ e (0)a x , (2.8) 

where the Hamiltonian Hjq(0) is also a modified JC 
one, given by replacing the parameters g and e in the 
modified JC Hamiltonian Hf c (e) with the scaled cou- 
pling constant g/e and the value respectively. While 
the Hamiltonians H-pc and Hf c (e) act in the state space 

T := C 2 ® L 2 (R), the modified JC Hamiltonian H% £ (0) 
acts in the chiral state space a x !F. Thus, we call the 
second term ea x Hf^ e (0)a x the chiral-counter Hamilto- 
nian of the modified JC Hamiltonian Hj C (e), and then, 

the pair {fff c (e) , ea x H% £ {0)a x } the chiral pair for the 
decomposition (|2.8p . We will show how to determine 
the parameterization for the modified JC Hamiltonian 
Hf c (e) as well as the naive meaning of the parameter e 
in Sec lIIIi We should note that Bermudez et al. have 
handled such a chirality to study a quantum phase tran- 
sition in the Dirac oscillator (4f| |4(| . We will find another 
phase transition in the chiral-counter Hamiltonian. 



We formulate our problem in the following. We denote 
each eigenstate of the modified JC Hamiltonian Hj c (e) 
by </3|(e), and its eigenenergy by E&(e), where the in- 
dex v runs over all integers: Hj G (e)<pf (e) = E&(e)ipf,(e), 
v = 0,±f,±2,---, with (^(e)|^(e)} = V Here let 
us take eigenvalues E&{e) so that the eigenvalue Eq(e) 
becomes the ground state energy when there is no inter- 
action (i.e., g = 0): E$(e) < E°(e) for v ^ 0. As shown 
in Sec lIIIi actually, we can take the index v so that the 
inequality, E^_, ,(e) < E^^(e), holds for every coupling 
constant g and the parameter e. As shown in Sec lIVl 
there is some chance that each of eigencnergies E&(e), 
v = — I, — 2, ■ • • , becomes the ground state energy when 
the coupling strength is large. Then, in the case where 
the eigenenergy E&{e) with the index v < becomes the 
ground state energy, it is called the superradiant ground 
state energy because of the physical mechanism found by 
Preparata |4l| - |44l ]. We denote the (superradiant) ground 
state energy of the modified JC Hamiltonians Hj C (e) and 

h jc(°) b y e jc( £ ) and E lc(°)> respectively. For these 
ground state energies, we define a quantity £i ow (g,£) by 

£ low (g, e) := Bf c (e) + E £$f(0). (2.9) 

Then, our main purpose is to approximate the ground 
state energy Epc of the Hamiltonian Hpc by using this 
quantity i?i ow (g,e:), and then, to study the role of the 
chiral part eE j^(0) in the ground state energy Epc- 
If the modified JC Hamiltonian Hf c (e) and its chiral- 
counter Hamiltonian ea x H^ (s)a x were commutable in 
the strong sense, i.e., [e ltH Jc( £ \ e *« r a>#jc , ( e )°'*] = for ev- 
ery i,s£M, the ground state energy Epc would be equal 
to the quantity E\ ow (g,s) which is then a constant func- 
tion of the variable e. But , unfortunately, they are not 
commutable in fact. Thus, we define the energy differ- 
ence £'diff(e) between the energies E-pc and E\ ovl {g,e) so 
that E FC = E iow (g,e) + E ditf (e): 

E diff (e) := E FC - E low {g,e). (2.10) 

Thus, the energy difference £'diff(£) represents the non- 
commutativity between the two component Hamiltonians 
of the chiral pair. Our first problem consists of estimating 
the energy difference i?diff(e) and computing the energy 

sEjg(0). The estimates of the energy difference £ , c iifi(e) 
clarifies how the non-commutativity shifts the quantity 
Ei ow (g,e) from the true ground state energy E F q. The 
estimate of the value of the energy eEf c £ (0) tells us how 
the effect of the chirality associated with the DP super- 
radiance is turned on and works as the coupling strength 
grows. 

For the normalized ground state ipE FC °f the Hamilto- 
nian Hpc we set transition probability amplitudes Af, (e) 
and Bl ,e {Q) as Af(e) := (^( £ )\^ Efc ) and Bg /e (0) := 
(<^ e (0) | t/'.Efc) respectively. Here the vector tf%/ £ (0) is 
the eigenstate of the modified JC Hamiltonian Hjg(0) 
with its eigenstate energy E^ e (0). Since we can prove 
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that the ortho-normal bases, {</?!(£)}> and {ff>^ e (0)}v, arc 
respectively complete in the state space T and the chi- 
ral state space o x T in the same way as proving Propo- 
sitions 3.1 and 3.2 of Ref. 4l|, we can expand the vec- 
tors ipE FC and a x ipE FC as ipE FC = Y. v A l{ e ) i P%{ E ) and 
a x ipE FC — J2 V B^ e (0)ipf,^ e (0). Inserting these expan- 
sions into the expression: Efc = (V'B PC |-ffpc|V'B F c} = 

(i>E FC \Hl c (s)\i()E FC ) + s(a x ipE FC \H^ e (0)\a x ->pE F c), we 
reach the following expansion: 



E FC = J2 \A*(e)\ 2 E*(e) + £ £ \B^(Q)\ 2 E^(Q). 

(2.11) 

We are interested in which term mainly contributes in 
this expansion. 

As precisely shown in Scc lIVl the ground state energies 
£j C (e) (resp. Ej^(0)) of the modified JC Hamiltonian 
becomes a superradiant ground state energy as the cou- 
pling constant g (resp. g/e) gets large. Thus, we indicate 

them as E^(e) := £f c (e) and £^(0) := E% E (0), re- 
spectively, with proper non-positive integers and i/**- 
Thus, to solve our problem of computing the ground 
state energy E^(0) (resp. Ej C (e)), we have to grasp 
the relation between the index i/** (resp. i>#) and the 
coupling constant g. By using the integers v„ and v*+, 
we have another expression of the quantity i?i ow (g, e) 

as Si ow (g,e) = (e) + eE^(0). Inserting Eq.^J^ 
into Eq. (|2.10j) with this expression, and combining the 
newly obtained expression of the energy difference and 
the equations, £„ \Af( £ )\ 2 = 1 and |B« /6 (0)| 2 = 1, 
we obtain the following expansion eventually: 

E dm (e) = J2 \Ai(e)\ 2 (EZ(e)-El(e)) 

+ e J2 \B^(0)\ 2 (E^(0)-E^(0)) 

(2.12) 

with El(e) > Esje) and E^ /e (0) > E^ £ r (0). Hence it 
immediately follows from this that the energy difference 
is non- negative: Edig(e) > 0. 

Define a lower bound £'ibd(e) and an upper bound 
S u bd(e) by -Bibd(e) := ei ow (g)-£ ; iow(g, e) and E uhd (e) := 
e upp(g) — -Eiow(g, s), respectively. Then, rewriting the in- 
equalities as ei ow (g) < E low (g, e)+E dm (e) = E FC < 
e U pp(g) leads to the estimates: 

max{0 , E lhd {e)} < E diS (e) < E uhd (e). (2.13) 

Namely, we obtain the inequalities: 

max{£'i ow (g,e) , E iow (g 7 e) + E lhd (e)} 
< E-pc 

<E low (g,e)+E uhd (e). (2.14) 

Then, our problem of estimating the energy difference 
E d is (e) means to tune the parameter e so that the upper 



bound E n \, d (e) becomes small for every coupling constant 
g. We can make our estimate better if we find a better 
upper bound than the bound e upp (g) and employ the 



-upper 



(g)- 



better one instead of the bound 

We make a remark on the order with respect to the 
coupling constant g in the estimates (|2.14[) : Actually, the 
order of the energy Ei ow (g, e) is given as Ei ow (g, e) ~ -g 2 
as g — ¥ oo, and then, we have the order 2 as: 



£i ow (g,e) +E diS (s) 



-g as g — > oo, 



(2.15) 



which is proved in SccjBj Thus, this is consistent with the 
order E F q ~ — g 2 (g oo) derived from the estimates 
(|2.4[) . Then, the following question arises: Each eigenen- 
ergy of the modified JC Hamiltonian has the order or 

1 (i.e., E$(e) g° and #| M (£) ~ ig 1 for v ^ 0. See 

Eqs. (|3~2]) and (pT3)) below). Why do the ground state 
energies E^ c {e) and Ej^(0), so the quantity Ei ow (g,e), 

have the order 2? i.e., £f c (e), E% £ (0), S low (g, e) - -g 2 ? 
We can find an answer for this question in Preparata's 
superradiance [H, HI, EE 13 • Mathematically, this or- 
der comes from the fact that the ground state energy of 
the modified JC Hamiltonian iJ| c (e) (resp. Hj^(0)) is 
made by the envelope consisting of all energy levels E& (e) 
(resp. £f /e (0)), v = 0,-1,-2,"-, as in Figs g] below. 
This envelope represents the appearance of the superra- 
diant ground state energies caused by the Dicke-type en- 
ergy level crossings. We use this fact in the ultra-strong 
coupling regime to approximate the ground state energy 
Efc by the quantity Ei OV! (g,e). On the other hand, we 
cannot use this fact in the weak coupling regime because 
the Dicke-type energy level crossings is special to the suf- 
ficiently large coupling strength. So, we employ scaling 
transformation to argue the order in the weak coupling 
regime later. 

We introduce the following notation. We call the pair 
of integers v* and i>** the superradiance indices, denoted 
by [ |^*| , ](e), v*, v** — 0, —1, —2, ■ ■ ■ . Namely, the 
superradiance indices \ \v* \ , \v^ \ \(e) mean the relations, 



E 



JC (e) = E^(e) and £f£ E (0) = E^'J(0). We note that 



the integer |z/*| (resp. |^**|) is the number of photons 
that the superradiant ground state is dressed with in the 
standard (resp. chiral) state space. Then, we will inves- 
tigate the effect of the chirality associated with the DP 
superradiance by finding the relation between the index 
f*» and the coupling constant g, and then, computing 

the the energy £-Ef£ £ (0) = sE^J(0). 

Our second problem is to estimate each transition 
probability |Af(e)| 2 from the ground state V'-Efc = 
Y, v Af,(e)(pf,(e) of the FC model to the state <pf( £ ) 
of the modified JC model, and also, each transition 
probability |£?^ e (0)| 2 from the ground state ^Bfc = 
E„ Bf /e (0)cr x ipl /£ (0) to t he state <r x <pl /e (0). Then, we 
can expect from Eq. (|2.12p that the transition probabil- 
ity |Ag(e)| 2 (resp. ji^^O)! 2 ) concentrates around the 



superradiant ground states: lim 



hm| 1 ,_iA„|- i .oo \B\ 



|f — | —¥oo 



\Al{e)Y 



or 



(0)| 2 = 0, which is argued in Sec lVIl 
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Thus, the number of photons which the ground state 
of the FC Hamiltonian is dressed with is dominated 
the states concentrating around the supcrradiant ground 
states. We will show at the tail end of Sec lVIl that the 
photon number is linked to the matrix size for the nu- 
merical computation of the ground state energy of the 
FC Hamiltonian. 



III. INTRODUCTION OF A KIND OF 
CHIRALITY 

To handle the counter-rotating terms in the ultra- 
strong coupling regime as well as in the weak or strong 
coupling one, we determine our parameterization of the 
frequencies w a and u) c . The naive meaning of the intro- 
duction of the parameter e is the following: Neither the 
rotating terms ftg(a^(T_ + cr + a) nor the counter-rotating 
terms hg(a)<j + + cr_a) can single-handedly make the sys- 
tem's energy spectrum. They need the free Hamiltonian 
to pay off. Thus, the parameter e represents how the 
rotating terms and the counter-rotating terms scrambles 
for the free Hamiltonian. Adopting the idea in Ref.(38j 
into our argument, we give the parameterized frequencies 
w a (e) and uj c (e) by 

LJ a (e) := (1 + e)cj a and cj c (e) := (1 - e)u c (3.1) 

for < e < 1. We employ this parameterization from 
now on. Then, we realize that the definition (|2.7[) of 
the modified JC Hamiltonian with the parameterization 
(|3.ip leads to the decomposition (|2.8j) . Here, we note that 
both energy spectra of the Hamiltonians ea x H^ e (0)a x 

and eHjg(0) are same, and that the equation, Hj C (0) = 
i/| c , holds. 

We now recall the eigenstates of the modified JC 
Hamiltonian H^ c {e) following the review in Ref . [l2| . 
Define states | \, n ) an d | \.,ri) for n = 0,1,2,- ■■ by 
| t,n) := |t)® \n) and | |,n) := \ i) <S> \n), respec- 
tively, for spin states, | \.) and | f), and the Fock states 
\n), n = 0, 1, 2, • • • . The state | f, 0) is an eigenstate of 
Hf c (e), i.e., 

J fff c ( £ )|t,0)=S t , o (e)|t,0) 
with its eigenenergy 

£ t ,o(e)=-^, (3.2) 

where A e :— (1 + e)uo a — (1 — e)uj c = A + e(cj a + cj c ) and 
Ao := cj a — lo c . We remember the assumption (|2.2[) . which 
implies the condition A e > for every e with < e < 1. 
For n = 0, 1, 2, • • • , set states |+, n)^ £ and |— ,n) g , £ by 

|+,n} giE := cos6»s(e)||,n) + sin6»s(e)|t,n + 1), 
\-,n) g , e :=-ame*(e)\l,n} + coae*(e)tf,n + l), 

where 6»|(e) := ifan" 1 (2g^T+T/A e ) if A e ^ 0, and 
6>«(e) = 7r/4 if A e = 0. Then, the states |+,n) gi£ and 



| — ,n) gy e are eigenstates of Hf c (e): 

Hl c (s)\±,n) e , s =El n (s)\±,n) s , e 
with individual eigenenergies 

El n (e) = (l-e)Hu c (n+l)±^A* + 4g2(n + 1) . (3.3) 

Therefore, we have the correspondence, Eq (e) = E^ (e), 
El H (e) = E'l n (e), and E s _ H (e) = E s _Je) with the re- 
lation \u\ = n + 1 for n = 0, 1, 2, • • • . The set of all ener- 
gies of fff c (e) is {£ t) o(e), E^ n {e) \ n = 0, 1, 2, • • • }. For 
simplicity, here, we denote the state |1\0) by |— ,— l) giE 
and | — , — l)g/ e ,o- Then, we recall the notation of the 
supcrradiancc indices: For integers := \ — 1 and 
n M := Iia^I — 1, the superradiance indices [|^*|, |z/**|](e) 
can be rewritten, and then, the equations, i?| c (e) = 

ElM) = E^ n ,{e) and E^{Q) = E^(0) = <^,(0), 
hold respectively. 

IV. SUPERRADIANT GROUND STATE 
ENERGY 

In this section, we recall the superradiant ground 
state energy of the modified JC model. The DP su- 
perradiance is a kind of phase transition so that the 
system of the 2-level system coupled to a photon-field 
switches from the initial ground state to a new non- 
perturbative one [H, [39|. According to the results in 
Refs.(4ll EU, there is the critical point of the coupling 
strength so that the following transition takes place: the 
energy E^_ _ 1 {s) = E^, (e) is the ground state energy 
of the modified JC Hamiltonian Hf c (e) for the cou- 
pling strength less than the critical point, but the energy 
E^L (s) replaces the old ground state energy E^_ _ 1 {e) 
and becomes the new one for the coupling strength more 
than the critical point. At the critical point, the modi- 
fied JC Hamiltonian has a degenerate ground state. The 
energies E*L n (e), n = 1, 2, • • • , also become the ground 
state energy in turn as the coupling constant g grows 
larger and larger, though they are primarily an excited 
state energy when the coupling constant g is very small. 
The new ground state energy E^_ n (e) is called the super- 
radiant ground state energy. The DP superradiance is 
a phenomenon special to very strong coupling strengths 
for the modified JC model as well as the standard JC 
model. The superradiant ground state energy E^_ „(e), 
ii = 0, 1, 2, • • • , cannot be obtained by the perturbation 
theory 

We can see the phase transition if we regard the ground 
state energy as a function of the coupling strength g: the 
energy E^ _ 1 (e) is a constant but the energy E^_ n {e), 
n = 1,2, is almost a first-degree polynomial func- 
tion as in Eqs. (|3.2[ ) and (|3.3[) . We can also explain 
this phase transition in terms of the photon number: 
the ground state is dressed with no photon in the case 
Ej C (e) = E*L _x(e)i on the other hand, it is dressed with 
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n + 1 photons in the case Ej C (e) = E^_ n (e) for each n 
with the condition n > 0. Because of the appearance of 
the supcrradiant ground state energy, there is surely an 
energy level crossing between the primary ground state 
energy and the superradiant ground state one, also, be- 
tween the previous superradiant ground state energy and 
the new one. This energy level cross ing is called the 
Dicke-type energy level crossing in Refs.[4l|,[42[ because it 
is due to the mathematical property special to the Dickc 
model (for instance, see the remark at the tail end of 
Sec.4.3 of Ref.jH). 

We can find such supcrradiant ground states and en- 
ergy level crossings for the standard and modified JC 
model by the numerical result in Figsf?] (a) and (b). 
Then, we can realize that the parameter e controls the 










(b) 





glh g/to of interaction (e=0) 



:nglh g/to of interaction (E-0.50) 



FIG. 4. Dicke-type energy level crossings among n {s), n — 
— 1, 0, 1, 2, • • - , of the (modified) JC Hamiltonian. Each color indicates 
individual index n of the energy n {c). Here cj a — uj c — oj. (a) 
standard JC model (e — 0); (b) modified JC model (e — 0.50). 

curvature of the ground-state-energy curve. Moreover, 
the ground state energy Ej C (e) (i.e., the envelope con- 
sisting of E^_ n (e), n = —1, 0, 1, 2, • • • ) makes the order: 
E^ c (e) ~ — g 2 as g — > oo as in Eq. (|2.15p . It has been 
shown in Ref s . |4l|, |42] how superradiant ground state en- 
ergies of the JC model appear in turn. So, we recall the 
rigorous description of this phenomenon here. 

To grasp superradiant ground state energy of our mod- 
ified JC Hamiltonians Hj C (e) and Hj{f(0) we introduce 
some constants. We can find a quantity G* (e) that com- 
pletely controls the crossing between the energy E^-,o(e) 
and the energy E^ n (e): 



ground state energy) Ej C (e) of the modified JC Hamil- 
tonian Hj C (e), we introduce a quantity G^(e) by 



2(1- £ )V (n + 1) 




for n = 0, 1,2, Then, we can state that for n = 
1, 2, • • • , the following holds: 

£ J s c ( £ ) = min{^, n _ 1 ( £ ),^.„( £ )} 

ifG+_ 1 (e)<g 2 <G+( £ ). (4.2) 

It is easy to check the inequality: 

G* n (e)<G+(e), n = 0,1,2,- ••, (4.3) 

when uj a = uj c . Then, the following question arises: 
which one is the supcrradiant ground state energy in 
Eq. (|4.2|) ? Before considering this question, we make a 
remark. Combining the crossing (|4.1j) and the superradi- 
ant ground state (|4.2p we have the following finding. In 
order that the energy E^_ n (s) becomes the superradiant 
ground state energy, it has to overtake not only E^.o(e) 
but also the energies E^_ m (e) with m < n. These are the 
Dicke-type energy level crossings. So, to answer to our 
question we have to consider the Dicke-type energy level 
crossing between the energies E^_ „_i(e) and E^_ n (e). 

In the case where u> = uj a = uj c and e — 0, it is easy 
to compare the two energies E^_ n _i(0) and E^_ n (0) as 
in Fig0ta) since A e = in this case: 



r^ lB _ 1 (o)<^, n (o) 

if and only if |g| < (V n + 1 + \/n)uj, 

< n _l(°) = <«(°) 

if and only if |g| = (yri+l + y/n)u>, 

£i,„_i(o) > ^,„(o) 

if and only if |g| > (yjn + 1 + y/n)uj, 



(4.4) 



G* n (e) := (l-ey<4 (n + 1) + 



A. 



1 



r)w c 



for n = 0, 1, 2, • • • . Namely, it gives a necessary and suffi- 
cient condition for the crossing between the two energies. 
Then, in the same way as in Ref. [4^, we can show the 
following facts: 



£ T ,o(e) < if and only if g 2 < G* n {e), 

E t>0 (s) = EZJe) if and only if g 2 = G* (e), (4.1) 
^ t ,o(e) > ££.„(e) if and only if g 2 > G* (e). 



Thus, the original ground state energy E^ : o(e) is not the 
ground state energy of the modified JC model any longer 
if the condition g 2 > G* (e) holds for each integer n. To 
seek the new ground state energy (i.e., the supcrradiant 



for 

Git. 



n = 0, 1,2, We here note the inequalities: 

! (0) = 4e 2 Lj 2 n < (VnTT + ^fn) 2 e 2 ' 



uj- < ie 2 uj 2 (n + 

1) = G+(0) for each n = 1,2,- ■•. Thus, the fol- 
lowing energy equations follow individually from the 
facts g2]) and (03]): £f£ e (0) 



= £f.^-i(0) provided 



e 2 G+_ 1 (0) = 4e 2 w 2 n < g 2 < (VrT+T + 0i) 2 £ 2 cj 2 ; and 



E% e (0) = E^ e n (0) provided (V^TT+ ^) 2 e 2 uj 2 < 
Ae 2 uo 2 {n + 1) = e 2 G+(0), for each n = 1,2, ■ • • . There- 



< 



fore, since e 2 G*(0) = e 2 u 2 = (VoTT + V0) 2 e 2 w 2 , we 
can conclude the energy equation giving the answer to 
the above question: 



EI(?(0) = E^(0) 



(4.5) 



for n = 0, 1,2,- ■ • if (\AT+T + y/n)' 
[Vn + ^+VnTlfe 2 ^ 2 - 



< 



< 
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It is not easy to compare the two energies n _ 1 (e) 
and n (e) under the condition e ^ 0. To see the cross- 
ing in such a case we introduce quantities G^ c (e) and 
Gl c (e) for n= 1,2,- •• by 

+ ^(" + 4 ) 2 + (^k) 2 } 

and 

respectively. Then, although we cannot give a necessary 
and sufficient condition for the crossing, we can give suffi- 
cient conditions for non-crossing and crossing. The same 
way as done in Ref.[42| brings us to the following facts: 

fg 2 < Gf(e) implies E^^e) < E*_Je), 
\g 2 > G s n c (e) implies > Eije). 

The first implication of the facts (|4.6p is an improvement 
of the result in Refs.flil |42| . Its proof is in SccJC] It is 
easy to check the inequalities: 

G™ c (e)<Gl c (e)<Gi(e), n = l,2,---. (4.7) 

Since we have the inequalities, G^ c 1 (e) < G^^ e (e), for 
each natural number £ with I > 2, if the condition g 2 < 
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G^F^e) holds, then we have the inequalities E^_ n (e) < 

E-,n+i( £ ) < E -,n+2^) < ■ < & tn +i(e) <■■ for each 
natural number I. On the other hand, since G^ c (e) > 
G^_ e (e) for each natural number I = 1, ■ ■ ■ , n — 1, if the 
condition G^ c (e) < g 2 holds, then we have the inequal- 
ities ££ >n ( E ) < E^ie) <■■■ < E^ n _ e (e) < ••• < 
E^ A (e). We note that Gl°(e) < G^(e). Therefore, we 
finally obtain the energy equation: 

£fcOO = E^ n (e) if G s n c (e) < g 2 < G^(e). (4.8) 

V. EFFECT CAUSED BY CHIRALITY AND 
NON-COMMUTATIVITY 

In this section let us set the condition, uj := u a = lu c , 
to make our argument simple, though we can argue our 
problem under the condition w a ^ u c with the assump- 
tion (|2.2p . Then, we prove that for every coupling con- 
stant g with < g < 1.193a;, there is a parameter e 
< e < 0.5 so that 

£iow(g,e) < E FC < E low (g,e) + QMhu, 
where the parameter e is determined when the coupling 
constant g is arbitrarily given. In addition, we show that 
there is a phase transition concerning chirality associated 
with the DP supcrradiance between the strong coupling 
regime and the ultra-strong coupling one: 

'eE% £ {0) = with v** = 

in the weak or strong coupling regime, 
eEjc(0) = huje\v^* \ — ?ig-\/|i / **| with z/** ^ 
in the ultra-strong coupling regime. 

More precisely, these statements can be proved as in Ta- 
ble HI where the function F(x,ri) of variables x and ?; is 
given by F{x, 7/) := |(1 — e~ 2x ) + n. 



coupling strength 


supcrradiance indices 
[\u.\, \v„\](e) 


chiral part 


upper bound of energy difference 
Bdiff(e) < Bubd(e) 


< g/w < E < 0.1 


[0,0](g/w) 





-Eubd(^) — F(g/uj, e)huj — hg 2 /uj with e — g/io 


0.1 < g/w < 0.2414 


[0. 1](0.10) 


o.whu - fig 


-Eubd(O.lO) = F(g/w, 0)fiw + fig - fig 2 /" < 0.25fiw 


0.2414 < g/w < 0.4 


[0, 1](0.20) 


0.20fiw - fig 


£ ubd (0.20) = F(g/u;,0)ftw + fig - fig 2 /" < 0.4fiw 


0.4 < g/ui < 0.5 


[0, 1](0.30) 


0.30fiw - fig 


F u bd(0.40) = F(g/w,0)ftw + fig - fig 2 /w < 0.5fiw 


0.5 < g/w < 0.9165 


[0, 1](0.40) 


0.40fiw - fig 


F ubd (0.40) = F(g/w, 0)fiw + fig - fig 2 /w < 0.53fiw 


0.9165 < g/u> < 0.9659 


[1,1](0.40) 


OAOhu - fig 


F ubd (0.40) = F(g/w, ^/0.16 + g 2 /w 2 )fiw + fig - Hg'/w < 0.53fiu; 


0.9659 < g/w < 1.193 


[1,1](0.50) 


0.50fiw - fig 


F ubd (0.50) < 0.56fiw 



TARLE I. Increase in Effect of Chirality and Estimate of Non-Commutativity of Chiral Pair 



We here remember the estimate (|2.14[) : Epc < 
E\ovt(g,£) + E u bd(s) = e upp (g). The phase transition as- 
sociated with the DP superradiance appears also for the 



modified JC Hamiltonian acting in the standard state 
space when the coupling strength approaches the deep- 
strong coupling regime, but Table fl] says that the phase 
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transition takes place faster in the chiral state space than 
it does in the standard state space. As for the coupling 
strength after g/uj = 1.193, although we have not proved 
yet, the numerical analysis in FigjSta) says the ground 
state energy Epc is bounded by the energy Ei ovr (g, 0.5) + 
0.56fiij from above. Also see FigJU The numerical anal- 
ysis in FiglS^b) says that the difference between the nu- 
merically computed ground state energy E of the FC 
Hamiltonian and the energy i?i ow (g, 0.5) + 0.50fouj is less 
than 0.025a;. From these numerical analyses as in Figs|5] 
the following problem arises: is the energy difference 
-Ediff (e) asymptotically equal to the vacuum fluctuation? 
We have not answered this question yet. 



1 0.0 





strength g/cn of inlerac 



strength g/G> of interaction 



FIG. 5. Set w a — uj c — d. (a) {£ — e upp (g))/?ia; (black solid line), 
where £ := -E low (g, e) + -Eubd(e) with e - g/hJ for < g/ui < 0.1; 
E = 0.1 for 0.1 < e/ui < 0.2414; e = 0.2 for 0.2414 < g/ui < 0.4; 
e = 0.3 for 0.4 < g/w < 0.5; e = 0.4 for 0.5 < g/u < 0.9659, and £ : = 
Bi ow (g, 0.5) + 0.56hu for 0.9659 < g/uj; (b) (£. ow (g, 0.5) + 0.50fi,a; - 
E)/huj (black solid line), where E is the numerically computed ground 
state energy. 



In the weak or strong coupling regime given by the 
region < g/u < e < 0.1, the supcrradiance in- 
dices [\v*\ , |i/**|](e) are [0, 0](e), and there is no effect 
of chirality associated with the DP superradiance, i.e., 
eEj^(0) = 0. Moreover, the ground state of the FC 
Hamiltonian is approximated as E-pc ~ E\ ow (g,e) = 
Ej C (e). Then, the ground state energy Epc is well 
approximated by the ground state energy Ejc of the 
standard JC Hamiltonian Hjq because of the limit 
lim e _j.o -Ejc ( e ) = -Etc • Here we note the following. The 
order of the energy E\ ow (g, e) with respect to the coupling 
constant g is given by -Ei ow (g,£) ~ — g° as g — > 0. Thus, 
the energy i?i ow (g,e) loses the dependence on the factor 
Hg 2 /uj c in the weak coupling regime, though the ground 
state energy Epc still keeps it as in Eq. (|2.3|) . We need a 
device. So, we introduce a scaling transformation (i.e., a 
renormalization) to argue this problem in Sec lVBI 

In the ultra-strong coupling regime, the index in 
the superradiance indices [\v*\ , |^**|](e) becomes non- 
zero (i.e., i>»» < 0), and the chiral part eE^(0) is turned 
on and works in E\ ovr (g, e). Thus, the effect of the chiral- 
ity associated with the DP supcrradiance appears. The 
effect of the chiral part increases more and more as the 
coupling strength grows large. In fact, Eqs. (|4.5[) and (|4.8[) 
say that as the coupling constant g grows large the in- 
dex increases faster than the index does. We 
find in Sec lVBI that for an arbitrarily fixed parameter 
e the index 1^*1 increases faster than the index in 
the superradiance indices , (See Tables HI1 

-Em. 



These argument say that a phase transition takes place 
concerning chirality associated with the DP superra- 
diance: the chiral part eEj^(0) is a constant in the 
weak or strong coupling regime, but it is almost a first- 
degree polynomial function of the coupling constant g in 
the ultra-string coupling regime. And moreover, in the 
ground state of the FC Hamiltonian Hpc , there is no pho- 
ton from the chiral-counter Hamiltonian sHjq (0) for the 
weak or strong coupling regime, on the other hand, some 
photons appear from the chiral-counter Hamiltonian for 
the ultra-strong coupling regime. In the case w a ^ w c , we 
can also find the phase transition concerning the chirality 
associated with the DP superradiance because the chiral 
part ei?j C e (0) is a constant, eE_ l x (0) = ~eh(u> li —uj c )/2, 
in the weak or strong coupling regime. On the other 
hand, it is almost a first-degree polynomial function, 
eE^! E n ^{Q) = ehu! c (n«« + 1) — Tig-y/n** + 1, in the ultra- 
strong coupling regime. 



A. Effect of chirality in i?i ow (g,e) 

Here we give some relations between the superradiance 
indices and the coupling strength. We pay our particular 
attention to how the index = n„ + 1 (i.e., how the 
effect of the chirality associated with the DP superradi- 
ance) increases as the coupling constant g grows for an 
arbitrarily fixed parameter e. To see that we consider suf- 
ficient conditions with respect to the coupling constant g 
for the index W^A. 



1. In Case [0,0] (e) 

We here consider the case where the supcrradiance 
indices are [0,0] (e). We note that there is no photon 
from the chiral-counter Hamiltonian in this case. The 
facts (|4.ip tell us that the necessary and sufficient con- 
dition for the superradiance indices [0, 0](e) is equivalent 
to the condition, < g 2 < G* (e) = (1 - e 2 )u 2 and 
< (g/e) 2 < Gq{0) = w 2 - This is also equivalent to the 
condition, 



0<g< 



if < e < 1/v^w 0.7071, 



Vl - e 2 uj ifl/ v / 2<£<l. 



(5.1) 



Therefore, we can obtain the weak and strong coupling 
regimes if we take the parameter e less than 0.1. 

When the integer n** — — 1 is equal to — 1, 

a straightforward computation shows that £i ow (g,£) = 

Ef c (e) (i.e., eE^ E (0) = 0) provided that w a = uj c . 
Therefore, there is no effect of the chirality associated 
with the DP superradiance in the weak or strong cou- 
pling regime. Actually, in the case u a ^= lo c , the chi- 
ral part depends on only the frequencies cj a and uj c : 
fs£j C e (0) = — eh(ij a — u c )/2. Thus, there is no effect 
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of the chirality associated with the DP supcrradiancc, 
either. 



2. In Case [0, v**](e) with v„ < 

We now consider the case where the supcrradiancc 
indices are [0,^**](e) with the condition i/** < 0. Wc 
note that some photons appear from the chiral-counter 
Hamiltonian in this case. We see here that the index 
+ 1 (i.e., the effect of the chirality associ- 
ated with the DP superradiance) increases as the cou- 
pling constant g grows. The facts (|4.1[) and (|4.5[) give us 
the sufficient condition for determining the superradiance 
indices. 

The sufficient condition for the superradiance indices 
[0, l](e) is implied by the condition: 

'(V2 + l)ew 

if < e < 72 - V2/2 ps 0.3827, 



the condition: 



euj < g < 



7T 



(5.2) 



if 72 - 72/2 < e < 1/72. 



Therefore, we can obtain the ultra-strong coupling regime 
if we take the parameter e more than or equal to 0.1. 

More generally, for each integer 77** > 1 (i.e., |z^** | > 2) 
we obtain the following fact. The sufficient condition for 
the superradiance indices [0,71** + l](e) is implied by the 
condition: 




< 



(7n** + 2 + 7n** + 
if < e < if (n** + 1), 

if K(n** +!)<£< K(n**), 



where 



K(n) := 



yj2(n + l) 2 - 2(n + 1) v^T + T) 



2(77+1) 



1,2,- 



These results say that we cannot expect the superra- 
diance indices [0,77** + l](e) under the condition 1 < 

77** and \/2 - 72/2 < e < 1 since K(n) < K(l) = 

for each natural number 77. On the other 
hand, we can expect the superradiance indices [0,n** + 
l](e) with 77** > 1 when parameter e are small enough: 
Let us arbitrarily take such a small parameter e that 



e 2 < 



1 



(V77** + 2 + v^7**+T) 2 + 1 



(5.3) 



Then, the energy equation (|4.5p says that a sufficient 
condition for the supcrradiancc indices [0,77** + is 



< g < (Vt7** + 2 + vn** + l)ew. 



(5.4) 



Meanwhile, for the supcrradiancc indices [1,1] (e) we 
obtain a sufficient condition for it by the Eqs. (|4.5p and 
KB as: 



'Vl~e^uj < g < (\/2 + i)ew 

if 0.38268 ps 72 - 72/2 < e < e , 
v/T - 7 2 "w<g<C 1 w C( e ) w 

if £ < £ < 1/72 0.7071, 



(5.5) 



2 w<g<C 1 wc (e)w 

if £o < e < £i, 

where Cj^ (e) is defined by 



Cf c (e) := ^(1 - e) {5(1 - e) + 741e 2 - 50e + 25} /2 , 

and the numbers Eq and ei with < So, £i < 1 ar e given 
by the solutions of the equations (72 + l)eo = Ci^^o) 
and £\ = C^ c (ei) respectively. The solutions e and £\ 
are respectively given as £o ps 0.496914 and £i ps 0.7500. 

As for the superradiance indices [1, 2](e), the Eqs. (|4.5| ) 
and (|4.8j) bring a sufficient condition for it: 



(5.6) 



(Vl~£^u < g < (7 3 + y/2)eu 
if 0.302905 ps 7l8 - 676/6 

< e < 72 - 72/2, 
(72 + < g < (73 + 72)ew 

if 72 - 72/2 < e < e 2 , 
(72 + l)ew<g<C2 vc (e)w 
if e 2 < e < e 3 , 



where the numbers e 2 and £3 with < £2, £3 < 1 are 
respectively given by the solutions of the equations (73+ 
72)e 2 = C] VG (e 2 ) and (72 + l)e§ = C™ c (e 3 ). Thus, the 
solutions are e 2 ps 0.424913 and £3 ps 0.496914. 

In the same way, we can seek a sufficient condition for 
other superradiance indices by using the results in Sec lIVI 



B. Estimate of Edis(e) 

Introducing the parameter e gives the degree of free- 
dom of the curvature to the energy-curve of the modi- 
fied JC Hamiltonian. For example, compare envelopes in 
Fig|4](a) and (b). In fact, wc obtain too many degrees of 
freedom of the curvature to determine the best parame- 
ter e. Thus, we determine here what is a better choice 
of the parameter e. To do that, we propose the criterion 
for the parameter e: when we fix the coupling constant 
g arbitrarily, we chose such a parameter s as the energy 
difference -Edifice) coming from the non-commutativity 
becomes as small as possible. 
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To determine the parameter e with this criterion, we 
employ the upper bound E u ^(e) defined in the estimates 
(|2. 13|) : -E'difi'(e) < £"ubd(£)- For simplicity we consider 
and examine only the candidates e = 0.10, 0.20, 0.30, 0.40 
and 0.5 in this paper. The upper bound -E u bd(e) is nu- 
merically estimated for some parameters e as in Fig|6] 
According to these numerical analyses, we should take 



2.5 : 

2 - 



o t- 1 1 1 1 1 

0.5 1 1.5 2 2.5 

strength g/co of interaction 

FIG. 6. Set cj a — uj c — ijj . i? u bd(e)/^- £ — 0.10(blue dashed-dotted 
line); e — 0.20(grecn dashed-dotted line); e — 0.30(green dashed line); 
e — 0.40(grocn dotted line); e — 0.50(bluc solid line); e — 0.60(magenta 
dotted line); e — 0.70(magcnta dashed line); e — 0.80(magcnta dashed- 
dotted line); e — 0.90(red dashed-dotted line). The envelope of all 
limes gives a better value of E u \,d(e) / hu) . 

the parameter e as e = 0.10 for < g/w < 0.2; e = 0.20 
for 0.2 < g/w < 0.4; e = 0.30 for 0.4 < g/oj < 0.5; 
e = 0.40 for 0.5 < g/w < 1.06; and e = 0.50 for 
1.06 < g/oj. We will prove this conjecture for the region 
< g/w < 1.19 (Tables HD— EIUl- The results suggests 
that it should be a good way that the rotating terms 
and the counter-rotating terms fairly share in the free 
Hamiltonian with each other. 



1. In Case [0,0] (e) 

In the case [0,0] (e): First, we can easily estimate the 
energy difference Edis(e) as < -Ediff(e) < (e + \ — 
ie _2g "/" 2 )fiw — Kg 2 /uj 2 in this case. We recall that the 
necessary and sufficient condition for the superradiance 
indices [0, 0](e) is equivalent to the condition (|5.1[) . Thus, 
by using this estimate together with Eqs. (|2.9[ ) and (|2.10l) . 
we obtain the approximation: 

E ¥C « E low {g,e), e<l, (5.7) 

because of the condition (|5.1j) . We remember that the 
weak and strong coupling regimes are obtained taking 
the parameter e less than 0.1 under the condition (|5.ip . 
and that the equation eEj(f(Q) = — eh{oj & — u) c )/2 holds. 
Thus, in particular, the ground state energy Epc has 
no effect of the chirality associated with the DP su- 
perradiance in these regimes provided that w a = lu c : 
E FC sa Ej C (e) by Eq. (|5.7|) . Moreover, we note that 
the ground state energy of the modified JC Hamilto- 
nian tends to that of the standard one as the param- 
eter e goes to 0: lira E _»o -^fc( e ) = -^fc because they 



are respectively given by Ej C = —h(uj a — lu c )/2 and 
Ej C (e) = Ej C — eh(uj a + w c )/2. Therefore, the ground 
state energy of the FC model is well approximated by 
that of the standard or modified JC model in the weak 
and strong coupling regimes. 

By the way, since the ground state energy of the mod- 
ified JC Hamiltonian Hf c (s) is a constant in the case 

= —1, it does not take care of the change in the cou- 
pling strength in the weak and strong coupling regimes. 
We try to cope with this problem in the following: We in- 
troduce a renormalized coupling constant g(e, a) for the 
parameter e and a new real parameter a as g =: e a g(e, a) 
to restore the order with respect to the coupling constant 
g in the weak or strong coupling regime. Forcibly adopt- 
ing the analogy of the standard QED in the weak cou- 
pling regime, the parameter e should be regarded as the 
rcnormalization constant and the parameter a should be 
1 /2. So, we set the renormalized coupling constant g W s(e) 
in the weak coupling regime as g W s(^) := g(e, 1/2). Al- 
though theory of QED for strong coupling regime has not 
been established yet, we employ g W s(e) of the standard 
strong coupling regime based on the estimates (|2.4p . 

Let us assume that we observe a coupling constant g ws 
with (g ws /w) 2 < l/y/2 for the ground state energy Epc 
in the weak and strong coupling regimes. We define an 
energy E^ a by E^ a := £^ c (g^ s /a; 2 ). Then, we obtain 
the estimates of the ground state energy E FC in these 
regimes as: 

- hglJuj = £fc < £fc < 0. (5.8) 

We will prove this fact in Scc lE II The estimate (|5.8[) 
says that the ground state energy Efc of the FC model 
can be well approximated with the renormalized ground 
state energy -Ef£ 3 of the modified JC model for the weak 
and strong coupling regimes (i.e., < g ws Cl). Thus, 
we should employ a sufficiently small parameter e = e ws 
in the weak and strong coupling regimes. Therefore, we 
obtain the following Table [TTJ 



superradiance indices 


region of coupling strength 


estimate of the ground state energy 


[o,o]( E ) 


< g/w < e < 1/V2 
(0 < g/w < 0.1) 


-l/y/2 < E fc /Hlo < 
(-0.1 < -gj,/w 2 = E^/hui < E FC /hu> < 0) 



TABLE II. Ground state energy estimates. 



2. In Case [0, i/„](e) with v„ < 

In the ultra-strong coupling regime, the chiral part 
eEjc(0) is turned on and works because the chirality 
associated with the DP superradiance appears. The ef- 
fect of the chiral part sEj^ (0) increases more and more 
as the coupling strength grows large because, as shown 
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in Sec IV A2l the index \v**\ of the superradiance indices 
increases faster than the index 

We estimate the energy difference Edis(e) in the case 
where n* = —1 and n** = 0, 1, 2, • • • . To do that, we 
only have to estimate the lower bound E\\, d (e) and the 
upper bound -E u bd(e) with the help of the inequalities 
(|2~T3l) and (|E2|) . 

We define two functions £ l °Yn ( x ) ancl £-i P n ( x ) 
by £-i, n ,A x ) ~ v 7 "** + 1 x -x 2 ' and ^^.(a:)* := 
\ + + 1 x — x 2 — \&~ 2x , respectively, for each in- 

teger n** = 0, 1, 2, ■ ■ • . Applying the mean value the- 
orem to the function f(x) = x 2 + \e~ 2x , there is a 
point xq between the non-negative points and x so 
that (/(0) - f(x))/(0 - x) = /'(.t ). Thus, we know 
that the expression | — x 2 — \ e ~ 2x is non-positive: 
5 — x 2 — \e~ 2x ^ = — 2cto(1 — e~ 2:r ") < 0, which implies 
the estimate £" pp „ ti ( x ) — V n ** -Fix. Let G-i,n,» be 
a positive constant given by a solution of the equation: 
2.t(1 — e~ 2x ) = y/n## + 1. We note that the positive 
solution is uniquely determined. 

Tra f/ie case [0, l](e): For the region (|5.2[) we have to 

assume that the parameter e satisfies e < l/y/2. We 
have the constant G-i,o ~ 0.745363 and the inequality 
£^ pp (G_i.o) < 0.53. Thus, the following estimates follow 
from the expressions (|E2j) : 

fc^o (I) 
< E dm (s) < faE^o (^) (5.9) 
< min{g/oj, 0.53}hu> —: £ Vm ,=\(e) 

in the region e < g/u> < mm{(y/2 + l)e, vl — e 2 } be- 
cause of the condition ()5.2|) . We note that £ , diff(e) < 
0.5hio for the coupling constant g satisfying the con- 
dition e < g/uj < min{(v / 2 + l)e, Vl — £ 2 , 0.6} since 
£- P i P o(x) < 0.5 for < x < 0.6. 

In the case [0,2](e): We put the restriction caused by 
the condition (|5.3[) on the parameter e: e 2 < {(\/3 + 
v^) 2 + l} -1 . We can compute the constant G_i i as 
G-\s w 0.88976 and estimate the value £ as 
S^ P A (G-i A ) < 0.87. Thus, the expressions (JE2J) come 
up with the following estimates: 

- ehu + fi^Ti (|) 
< E diS {e) 

< -ehu + hwE^ (i) (5.10) 

< (min{\/2g/u; , 0.87} - e^j hw 

=:^..=a(e) 

in the region (-\/2 + l)e < g/cj < (V3 + V2)e, which is 
determined by the condition (|5.4p . 



Tra the case [0, 3](e): We put the restriction caused by 
the condition (|5.3j) on the parameter e: e 2 < {(2 + 
•\/3) 2 + I} -1 - The constant G_i.2 is obtained as G_i.2 = 
1.000964 and the value £" pp 2 (G_i2) is estimated as 
f" pp 2 (G_i :2 ) < 1.18. Thus, we can derive the following 
estimates from the expressions (|E2[) : 

- 2ehw + hw£ x °l 2 (i) 
< E dm (e) 

< -2ehu + fiw£"P p 2 (i) (5.11) 

< (min{V3g/a; , 1.18} - 2ej hw 

=: £ v ..=3(e) 

in the region (\/3 + \/2)£ < g/w < (2 + \/3)e which is 
implied by the condition (|5.4[) . 

in i/ie case [0,4](e): We put the restriction caused by 

the condition (|5.3[) on the parameter e: e 2 < {(v5 + 
2) 2 + 1} _1 . The constant G_i. 3 is computed as G_i :3 = 
1.09837 and the value £ "^(G-i^) satisfies the inequal- 
ity £^ p 3 (G_i, 3 ) < 1.47. Thus, by the expressions (|E2]t . 
we obtain the following estimate: 

- 3efiw + ftw£i7 3 
< ^diff(e) 

< -Sehoj + hw£ ^ pp 3 (i) (5.12) 

< (min{2g/w, 1.47} - 3e)fe»; 

=: £ v „ = 4:{e) 

in the region (2 + \/3)e < g/w < (V5 + 2)e coming from 
the condition (|5 ,4[) . 

Next we estimate the energy difference E d is(e) in the 
case where n*, n** = 0, 1, 2, • • • with the help of the in- 
equalities (pTT3"l) and (|E4|) . 

Define two functions £i° w „„(x) and ^ pp n „(x) by 

4° w n»» •= vk + i) x 2 + £ 2 - (n. + 1) + f!_°r;„„(x) 

and ^ pp „(x) := ^n, +l)i 2 + £ 2 - (n, + 1) + 
^•-i P n»»( a; ) respectively. Then, we note f^Jj (x) < 
v /K + l)x 2 + £ 2 - K + 1) + ^ P 1 P „^(G_ 1 ,„,J** 

In the case [1,1] (e): We put the restriction on e as 

0.38268 w V2 _ V2/2 < e < e x w 0.7500 as in the re- 
gion (|5.5j) of the coupling strength. Since there is a point 
Go.o so that £qq P (Go,o) = max x £q q P (x), the estimates 
£12} lead to 

(|) < ^diff(£) < fa)S%ff (^) (5.13) 
< ^£ upp (G ,o). 
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Wc note here that £o£ p (G ,o) < 0.53 for e = 0.40 and 
£o P p {G ,o) < 0.56 for s = 0.50. 

In the case where the coupling constant g is about in 
the region 0.1 < U) < 1, we obtain Tables Hill - fVl by the 
estimates ()5.9|) - (|5. 13[) . Then, the index is just but 
the index z/** grows. 



supcrradiancc indices 


region of coupling strength (e — 0.10) 


estimate of energy difference 


[0,1](0.10) 


0.1 = e < g/u < (V2 + 1)£ S3 0.2414 


£ diff (0.10)//iw < g/w < 0.2414 


[0,2] (0.10) 


0.2414 S3 (V2 + l)e 

< g/w < (V3 + V2)e S3 0.3146 


£ d iff(0.10)/ftw < V2g/w - 0.1 < 0.3449 


[0,3] (0.10) 


0.3146 S3 (V3 + V2)e 

< g/w < (2 + \/3)e S3 0.3732 


£ d iff(0.10)/ftw < V3g/" - 0.2 < 0.4464 


[0,4](0.10) 


0.3732 S3 (2 + V3)e 

< g/w < (VE + 2)e S3 0.4236 


E di{{ {0.W)/hu < 2g/w - 0.3 < 0.5472 


TABLE III. Energy difference estimates {e — 0.10). 


supcrradiancc indices 


region of coupling strength (e — 0.20) 


estimate of energy difference 


[0, 1](0.20) 


0.2 = e < g/w < (V2 + l)e S3 0.482843 
(0.2414 < g/w < 0.4) 


£ d iff(0.20)//iw < g/w < 0.482843 
(Bdiff(0.20)/fto) < g/w < 0.4) 


[0, 2] (0.20) 


0.482843 S3 (V2 + l)e 

< g/w < (V3 + \/2)e S3 0.62925 


£ d iff(0.20)//iw < £„„»=2(0.20) < 0.67 


[0, 3](0.20) 


0.62925 S3 (V3 + V2)s 

< g/w < (2 + V3)e S3 0.74641 


£ d iff(0.20)//iw < £„„»= 3 (0.20) < 0.78 


[0,4] (0.20) 


0.74641 S3 (2 + V3)e 

< g/w < (V5 + 2)e S3 0.84721 


£ d iff(0.20)//iw < £„„ =1 (0.20) < 0.87 


TABLE IV. Energy difference estimates (e — 0.20). 


supcrradiancc indices 


region of coupling strength (e — 0.30) 


estimate of energy difference 


[0, 1](0.30) 


0-3 = e < g/u: < (v / 2 + l)s S3 0.724264 
(0.4 = e < g/w < 0.5) 


£ d iff(0.30)//iw < £„,» = i(0.30) < 0.53 
(-E d iff(0.30)/?iw < g/u, < 0.5) 


[0, 2](0.30) 


0.724264 S3 (V2 + l)e 

< g/w < (V3 + \/2)e S3 0.9438793 


,E d iff(0.30)//iw < £„„»=2(0.30) < 0.57 



TABLE V. Energy differ ence estimates (e — 0.30). 

On the other hand, in the case where the coupling 
constant g is more than about 10, we obtain Tables IVII 
and I VIII by the estimates (|5.13p . Then, the index is 
beginning to be 1. 
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supcrradiancc indices 


region of coupling strength (e — 0.40) 


estimate of energy difference 


[0, 1](0.40) 


0-4 = e < g/w < VI - e 2 S3 0.9165 
(0.5 < g/w < 0.9165) 


,E d iff(0.40)/riw < = ! (0.40) < 0.53 
(B d iff(0.40)/riw < £„„„ = i(0.40) < 0.53) 


[1, 1R0.40) 


0.9165 S3 VI - e 2 < g/ui < (V2 + l)e S3 0.9659 


B d iff(0.40)/fiw < 0.53 



TABLE VI. Energy difference estimates (e — 0.40). 



supcrradiancc indices 


region of coupling strength (e — 0.50) 


estimate of energy difference 


[0, 1](0.50) 


0.5 = e < g/w < VI - e 2 S3 0.866 


£diff(0.50)/fiw < £„,,„ = i(0.50) < 0.53 


[1,1](0.50) 


0.86603 S3 VI - £ 2 < g/w < Ci(0.5) S3 1.1938972 
(0.9659 < g/w < 1.193) 


£diff(0.50)/frw < 0.56 



TABLE VII. Energy difference estimates (e — 0.50). 



Following the numerical analyses in Fig|6j the case 
where e = 0.40 gives a good energy-difference estimate 
for the region 0.5 < g/us < 1.06 among the candidate 
e = 0.10,0.20.0.30,0.40, and 0.50. But, as in Table ED 
our bound obtained is 0.9659. It is caused by the fact 
that the bound C^ vc (e) is not best possible. 

Using the results in Tables E]- IVlTl we can conclude 
the results on the estimates of the energy difference in 
Table H 



VI. CONCENTRATION OF SUPERRADIANT 
GROUND STATE 

We estimates the transition probabilities |^4s( e )| 2 anc j 
|-Bf(0)| 2 under the condition ui = w a = cu c in this sec- 
tion. For simplicity, we reset the notations of the tran- 
sition probability amplitudes A^{±,e) and £?^ 6 (±,0) 
as Az(±,e) := £ , g (n, ±|^ FC ) g , £ and B?/ £ (±,0) := 
o, g / £ (n, ±\a x ipE FC )g/efi, respectively, for n = 0,1,2,- ••. 
Here we denoted the ground state of the FC model 
by the symbol ipE FC - We meant the state | t, 0) by 
|— , — l)g, E and | — , — l) g / £! o- So, for the state we set tran- 
sition probability amplitudes A^ 1 (—,s) and B^(—,0) 
as Ai^-.e) := e , g (-l, -|^ FC ) g , £ and £^(-,0) := 
o,g/ £ (-l J -k:i : ^Epc)g/E,o, respectively. 



A. Weak and Strong Coupling Regimes 

We will show the following estimates in SeclDl 

1 - 4 < l^xC- e)| a + ^(-^O)! 2 < 1. (6.1) 
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Since the FC Hamiltonian Hpc converges to the free 
Hamiltonian (huj/2)a z + hg{a*a + 1/2) as g — > in the 
norm resolvent sense, combining Lemma 4.9 of Ref.(6l| 
or Theorem VIII.23 of Ref.[62| with the fact (|ETJ). we 
have the limits: 



g->0 



1 and 



lim \B' 

g-s-o 



i(-,o)r 







in the weak coupling regime determined by the condition 
(|El|) . Therefore, the estimates (|6~Tj) and the limits 



show the limit: lim 



g^o We 



- |t;0). Namely, we have 
proved that the ground state of the FC model is well 
approximated with that of the (modified) JC model in the 
weak and strong coupling regimes as well as the ground 
state energy. 



B. Ultra-Strong Coupling Regime 



Eq. (|2.12[ ) leads to the inequalities: 



0<K(±, £ )|^<min 1, 



0< |BS/ £ (±,0)| 2 < min 1 



E diS {e) 



e(E s J E n (0) 



<L,(0)). 



(6.3) 

Let us fix any superradiance indices [ ](e) = 

[n„ + Ln„ + l](e) now. Then, we have the limit 



E - El c (e) | = Iim ]n _ n . Hoo \E - E* t (e) \ 



Iim| n _ n . 1^.00 

oo for energy E = E± n {e) by using the concrete ex- 
pressions of the eigenstate energies E^ :0 (e) and E± n (e). 

Similarly, we have the limit Hmi n _,„ li j_ >00 \E—E^q 1 ' 



limi, 



-E%(l(0)\ = oo for energy E = £± 7 * (0). 



A |n— ti**|— >oc 

Since the energy difference E d is(e) is independent of 
the natural number n, we can conclude from inequal- 
ities (|6.3[) that the transition probabilities |>lj> (±, e)| 2 
and |_B^ e (±,0)| 2 concentrate around the superradiant 
ground states: 

lim K(±, £ )| 2 =0, 

|n— n* \—>oo 

lim |5^(±,0)| 2 =0. 

_ |n— n** 1— foo 

Therefore, the number of photons that the ground state 
of the FC model is dressed with is determined by the 
number of photons that the states around the superradi- 
ant ground state of the modified JC model are. 

At the tail end of this section we give a remark on 
the accuracy of the numerical analysis of the ground 
state energy of the FC model. The fact (|4.2p says that 
the superradiant ground state energy of the modified JC 
Hamiltonian Hf c (e) appears for the coup ling constant 
g \/ ti* as Preparata predicted in Ref.|43j (also see 
Ref.Jij]). Combining inequalities p.!4[) and (|4~2j) , we 
can say that the superradiant ground state is the state of 
the quasi-particle dressed with about photons when 
the coupling strength is so strong that the inequalities, 



G^_ 1 (e) < g 2 < G+ t (e), hold. Let us fix any coupling 
constant g satisfying these inequalities, and denote it by 
g„ t . Namely, G^__ 1 (e) < g 2 , < G+^e), which implies 
the equation: 



lim — : 



4(1 -e)W c 



for the modified JC Hamiltonian Hf c (e) 
have the limit 



(6.4) 



Similarly, we 



t2 



lim 



(6.5) 



for the modified JC Hamiltonian Hj^(0). Let us take 
0.5 as the parameter e and assume the condition w := 
Lj a = lo c . Then, the concentration of the superradiant 
ground state together with the limits (|6.4p and (|6.5[) says 
that for the given coupling constant g we have to han- 
dle an n x n- matrix with n ss g 2 /w 2 when we numeri- 
cally calculate the ground state energy of the FC model. 
Namely, if we use the n x n-matrix for the numerically- 
calculated ground state energy, the energy becomes in- 
accurate after the about dimensionlcss coupling strength 
g/oj = y/n. For example, see FigsJT] (a) and (b). As in 




FIG. 7. Set 



3 (g)//iu(rcd solid line), numerically- 



calculated ground state energy (blue dotted line), (a) n — 100 (\/n — 
10); (b) n = 200 (Vn S3 14). 

the inequalities (|2.4[) the ground state energy Epc has to 
be less than or equal to the upper bound e upp (g). How- 
ever, the numerically-calculated ground state energy has 
inaccuracy like e upp (g) < Epc for the range of the cou- 
pling strength after the dimensionless coupling strength 
g/oj = jn. 



VII. CONCLUSION 

We have proposed a Hamiltonian formalism by de- 
composing the FC Hamiltonian into the modified JC 
Hamiltonian and its chiral-counter Hamiltonian. The two 
Hamiltonians play the chirality in the FC Hamiltonian. 
We have characterized the ground state energy of the FC 
Hamiltonian in all regimes by using our decomposition 
and estimating the non-commutativity between the two 
decomposed Hamiltonians. That is, we have shown that, 
in the weak or strong coupling regime, the ground state 
energy of the FC model can be well approximated by 
that of the (modified) JC Hamiltonian. In this process, 
there is no effect of the chirality associated with the DP 
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supcrradiance. On the other hand, we have shown that, 
in the ultra-strong coupling regime, the effect of the chi- 
rality is turned on and works, and moreover, it increases 
as the coupling strength grows large. In this process, 
a phase transition concerning the chirality mathemati- 
cally appears when the atom-cavity system plunges into 
the ultra-strong coupling regime from the strong coupling 
one. If the phase transition is experimentally observed, it 
means that the DP supcrradiance is indirectly observed. 

We have employed the upper bound e upp (g) for the 
ground state energy Epc of the FC Hamiltonian in this 
paper to estimate the energy difference £'diff(e) between 
the energies E F c and Si ow (g, e). The latter is given by 
the sum of the ground state energies of the modified JC 
Hamiltonian and its chiral-counter Hamiltonian. In our 
method, however, if we can find a better upper bound 
than the bound e upp (g), we can make the estimate better 
by using the better bound instead of the bound e upp (g). 
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Appendix A: Proof of Eqs. (|2T5|) and (|276 



The inequalities (|2.4p immediately lead to the 
limits for the asymptotic behavior of the ground 
state E-pc'- hrn g ->.o -Efc = -fi(w a — lu c )/2 and 
lim g _ i . 00 i?Fc/(% 2 / w c) = 1- The former and latter re- 
spectively imply rhe first and second of Eqs. (|2.5p . 

To prove the third of Eqs. (|2.5j) and E g (12.61) . we follow 
the method to prove Theorem 10 of Ref.[55j on the spin- 
boson Hamiltonian. As in Eq.(2.14) of Ref . [Ho} , using 
the unitary operator Uq := (l/\/2)(co — ify), we have 
the representation: 



U„H FC U = hhj c (^a + + hg(T z (a + a + ) 



We apply the perturbation theory to the operator 
U^HpcUo/Hg. Namely, we regard the term — (u a /2g)cr x 
as a perturbation for the operator (uj c / g)(a) a + 1/2) + 
a z (a + a + ). Then, the infimum eigenvalue uj c /2g — 
l/(w c /g) splits into two eigenvalues: w c /2g— l/(w c /g) ± 
(o; a /2g)e- 2/( " c/s)2 + o(ul/g 2 ). Multiplying this by % 
and judging from the inequalities (|2.4[) . we obtain the 



expressions: 



E 





hiO c 


fig 


FC — 


~ ~ 








UJ C 


1st 


fllij c 


ng 


FC — 


~2~ ~~ 





hui a 
~Y 

~2 



± e -2(g7w c ) 



Therefore, we eventually reach the third of Eqs. (|2.5[ ) and 
Eq. (|2.6p . The former says that the limit G(g) — > holds 
as g — > oo. 



Appendix B: Proof of Inequality (I2.15[) 

Applying the inequalities (2.66) in Lemma 2.2 of 
Ref. [38l] . we have the following order estimates: 



J_ <i imm f|M2 

— £ g^oo hg Z /UJ c 



EUe) 
< limsup . J P V ; < 



1 • c £E JC 

- < hminf - ■ , . 

e g^oo ng^/u}, 



hg 2 /uj c 
/£ (0) 



4(1 -s)' 



(Bl) 



eE% £ (0) 
< limsup - J 9 L . 
g^oo ng z /uj c 



< . 

~ As 



These estimates imply the following inequality: 

71 T - hmmt fe 2 / 

e(l - e) g->oo ng^/ojc 



, ,. ^low(g,£) . 
< limsup < 

g^oo ng z /U} c 



1 



4 £ (l-e) : 



(B2) 



Meanwhile, combining the definition of the upper 
bound e upp (g) in the inequality (|2.4[) with the inequal- 
ities (|2.13p and (jBip . we reach the following bound: 



< limsup < 



g->oo 



hg 2 /u 



e(l -e) 



(B3) 



Thus, we know that the order of the energy difference 
Sdift(e) is given by E diS (e) = 0(g 2 ) as g -> oo. T here- 
fore, we eventually obtain the order estimate (|2.15l) . 



Appendix C: Proof of the 1st of the facts (14. (i 



We give a proof of the first statement of the facts (|4.6j) 
here. Set the number 5 as S := (1 — e)uj c for simplicity. 
Define a function g(r) of a real variable r by g(r) := r 2 — 
(4n+l)r— (A £ /5) 2 . Then, the two solutions r± of g(r) = 



0arer ± = (2n + 1) ± (2n + |) + (A E /S) 2 . Since one of 

them is negative (i.e., r_ < 0), if r satisfies the condition 
< r < r + . then g[r) < g(r + ) = 0. Inserting (g/5) 2 
into r in the above, we have the following statement: the 



condition (g/,5) 2 < 2n + \ + ^(2n + \) 2 + (A e /,5) 2 im- 
plies the inequality (g/5) 4 < (An + l)(g/<5) 2 + (A E /S) 2 . 
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Multiplying both sides of the second inequality in the 
above by 4(5/g) 2 , the second one implies 4(g/<5) 2 < 
4(4n + 1) + 4(A £ /g) 2 . The right hand side of this is 
decomposed and bounded as: (A £ /g) 2 + 4(n + 1) + 
(A £ /g) 2 + 4n + 2{( A £ /g) 2 + 4n} < (A E /g) 2 + 4(n + 1) + 

(A £ /g) 2 + An + 2 V(A £ /g) 2 + 4(7i + 1) V(A E /g) 2 + 4n = 
{ v /(A £ /g) 2 + 4(n+l) + ^(A £ /g) 2 + An } . Eventu- 
ally, we obtain the mathematical statement: 

(f) 2 < 2 » 1 




Meanwhile, it is easy to show that the mathematical 
statement: 

^-.n-i( £ ) it £"-.n( £ ) i s equivalent to 




where the symbol jj denotes the mathematical symbol 
<,=, or >. Combining the two statements (|C2[) and 
(|C2[) , we can complete the proof of our desired result. 



P Q >1-N. 

By Lemma 4.3 of Ref.[6l[, we have the inequality: 



< 



(Dl) 



(D2) 



where the vector V'Efc was the normalized ground state 
of the FC model. It is proved in the following. Us- 
ing the commutator [i/pd a]^£ F c = {Hfc — £ , Fc)aV , E F c 
with [H-pQi o] = ^huja — hga x , we reach the so-called 
pull-through formula: aipE FC — ~ %(#fc — -^fc + 
huj)~ 1 <T x 'il)E FC ■ Applying this pull-through formula to the 
term (tl>E FC l a HV'-Efc) and using the Schwarz inequality, 



we obtain the inequality (|D2[) . The two inequalities (|D1 
and (|D2[) imply the lower bound: 



(V ; B F d cr O-Po|V'£;Fc) > 1 



(D3) 



On the other hand, the Schwarz inequality brings the 
upper bound: 



(lpE FC WoPo\^E FC ) 

< (iPefcHefc) 1 ^ 2 (<7oPo4>e fc \°'oPo4> Efc / 

= {^E FC \lpE FC ) = 1. 



(D4) 



Since the equation lA^-,^ 2 + \B S J^{-, 0)| 2 = 
("0-Efc l fJ o-Po|'/'-B F c) follows from the straightforward com- 
putation, we eventually obtain our desired estimate (|6.1j) 
by the bounds (|D3)) and JD4|. 



Appendix E: Proof of Estimates _Ediff(e) 
1. In the Case [0,0] (e) 



In the condition (|5.1[) equivalent to the superradiance 
indices [0, 0](e), the following estimates are derived from 
the inequalities ([23]), (f2TT3|) . and (|2~T4j) : 



£;f c (£)=Ei ow (g ) £)< J E PC <0 J 
< g < ew with < e < l/\/2. 



(El) 



Since the equations -Eiowfee) = Ej C (e) = ~ £ ^ L0 hold 
for [0, 0](e), it is immediately follows form the definition 
of the upper bound -E u bd(e) and the estimates (|2.13[) that 



Appendix D: Proof of estimate (|6.ip 

In this appendix we give a mathematical proof of 
the estimate (|6.1[) from the below. We define the or- 
thogonal projection operator P n by P n ip := (n\^)\n) 
for the normalized photon state |n). Then, we have 
the expression of the number operator N := a^a as g 
N = Y^=n n Pn = Y^=i n Pn, so that the equation 
Po + N = P + nPn, which implies the operator 

inequality Pq + N > Y^=o P< 1 = 1- Thus, we reach the 
operator inequality: 



< Sdifi(e) < E uhA (e) = (e+|- 



7« 



The region of the bare coupling constant g in the es- 
timates (|E1[) leads to that of the renormalizcd coupling 
constant g ws (e) as < gws(e) < w ^h < e < 1/V%- 
Since we assumed that the observed coupling constant 
g ws satisfies the condition (g ws /ui) 2 < 1/V2, the quan- 
tity e ws satisfies the condition < e ws < l/\/2- Since 
?w S (e ws ) = gws, we have the equation g = ^/s^g ws = 
g^ s /tu. Meanwhile, we have the equation e ws cu = g^ s /w. 
Thus, the quantity e ws satisfies the last condition of the 
fact (|E1|) , and we can take e ws as the parameter e, i.e., 
£ = £ws- Then, the first part of the fact (|E1[) says that 
—hg^/Lj = -e ws fruj = i?| c (e ws ) < E FC , which is the 
lower bound of the estimates (15.81). 



2. In the Case [0, v««](e) with v„ < 

In this case we have the value of integers n* and 
as: = —1 and n*, = 0,1,2,- ■-. Using the concrete 
expressions of the lower bound -Eibd(e) and the upper 
bound -Eubd(e), we have their other expressions: 



E\ hd {e) = hu -en** + ~Vn** + 1 - ^ 

L U! U! z 
r 1 g 

E u bd(s) — fru 77 - H V"** + 1 - 

L2 cj 

_1 -2g 2 /^ 



T; n* 



0,1,2, 
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with the help of the inequalities (|2.13[) . Then, these equa- 
tions lead to the inequalities: 



-e (n** + 1) 



fko. 



(E2) 



The inequalities ([E2j) imply the results, (|5l?|) . (|5.10p . 
(I5TT1) . and (I5T21) . 



In the case where n^n** = 0, 1,2, similarly to 
the above, we obtain the expressions of the lower bound 



-Eibd(e) and the upper bound i? u bd(e) as: 
-Eibd(e) = fiw — (n* + 1) - e(n*, - n*] 



£ 2 + ^(n* + l) 

2 



-E ; ubd(e) = ^ i — (n* + 1) - e(n** - n*) 



£ 2 + ^(n, + l) 



+— + 1 



(E3) 



Then, using the functions and Eqs. (|E3[ ). we can estimate 
the lower bound Ei\,d(s) and the upper bound E u ^(e) 
as: 

'e«w(n. - n**) + 7^47n«, f 1 ) < ^ibd(e), 

' W /K\ ( E4 ) 

^ubd(e) < eftw(n, - n..) + ^^ p n „ (^J 

for the superradiancc indices [ \v*\, \v**\ ](e) = [n* + 
l,n*» + l](e), which implies the result (|5 . 13|) . 
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